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Marangoni propulsion is a form of locomotion wherein an asymmetric release of surfactant by 
a body located at the surface of a liquid leads to its directed motion. We present in this paper 
a mathematical model for Marangoni propulsion in the viscous regime. We consider the case of a 
thin rigid circular disk placed at the surface of a viscous fluid and whose perimeter has a prescribed 
concentration of an insoluble surfactant, to which the rest of its surface is impenetrable. Assuming 
a linearized equation of state between surface tension and surfactant concentration, we derive ana- 
lytically the surfactant, velocity and pressure flelds in the asymptotic limit of low Capillary, Peclet 
and Reynolds numbers. We then exploit these results to calculate the Marangoni propulsion speed 
of the disk. Neglecting the stress contribution from Marangoni flows is seen to over-predict the 
propulsion speed by 50%. 



I. INTRODUCTION 

a: 

. The study of animal locomotion, as carried out by zoologists and organismal biologists [H, [Ij , has long been a source 
of new problems in fluid dynamics, either because fluid flows pose biological constraints which deserve to be quantified 
Q or because biological situations lead to interesting and new fluid physics [1, Hj. Broadly speaking, there exist two 
types of self-propelled motion in fluids. In the flrst type, locomotion is induced by shape changes. This category 
includes all classical biologically-relevant work at both high and low Reynolds numbers, and encompasses flying birds 
al and insects 16l-l9i, swimming fish aEp, the locomotion of microorganisms jl2| - |l5{ and interfacial propulsion 
^ ■ M- 

The second class of self-propulsion is one where motion occurs without relative movement between the body and the 
fluid. In that case, transport arises due to specific interactions between the body and its surrounding environment. 
Beyond the trivial case of force-driven motion such as sedimentation (iTj . this category includes transport driven 
by gradients, either externally applied [l8l - l2l| or self-generated [22| - [26| . In this paper, wc consider a specific setup 
belonging to this second category, namely locomotion at the surface of a liquid driven by self-generated surface tension 
gradients, a situation sometimes referred to as Marangoni propulsion [l6| . 

A gradient in surface tension at the surface of a liquid can give rise to transport for two distinct physical reasons. 
First, surface tension acts as a force along contact lines. A body subject to a surface tension gradient with nonzero 
average will thus experience a net force, and will have to move. Well-studied examples include droplets moving 
on substrates with prescribed gradients in surface energy [l^, HH, H^] , or reacting with a substrate to self-generate 
similar asymmetries [2^, H^. Surface tension gradients can also lead to transport because of Marangoni flows. The 
tangential stress condition at the interface between two fluids states that tangential viscous stresses arc in mechanical 
equilibrium with gradients in surface tension, and thus any imbalance in surface tension leads to a flow (30j . Examples 
of Marangoni transport include the motion of droplets and bubbles in thermal gradients [l8l - l20j and droplets subject 
to internal releases of surfactant [3l| . 

Marangoni propulsion refers to the situation where a body, located at the surface of a fluid, generates an asymmetric 
distribution of surface active materials, thereby prompting a surface tension imbalance and a Marangoni flow, both of 
which lead to locomotion. A simple illustration of this phenomenon is illustrated in figure [T] using a toothpick dipped 
in a low surface tension dishwasher liquid. The surface tension at the contact line near the rear of the toothpick 
is smaller than that near the front, which leads to forward propulsion. Marangoni propulsion in Nature, recently 
reviewed by [l^, is of relevance to the behavior of insects. The most studied is Stenus (rove beetle), which is observed 
to recover from a fall on a water surface by quickly releasing surface active materials j33 - l3a | . Marangoni locomotion 
is also used by Dianous (rove beetle) [s^l and Velia (small water strider) [s^l • 

Experimentally, a lot of work has been carried out on synthetic Marangoni propulsion, mostly regarding the dy- 
namics of so-called Camphor boats, where the dissolution of Camphor grains at the surface of water leads to their 
propulsion [bsI . [39j . For that system, the Marangoni flow has been visualized [40j , and one-dimcnsional modeling 
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FIG. 1: Example of Marangoni propulsion in the kitchen sink. The tip of a toothpick is dipped in dishwasher liquid, and is 
then deposited at the free surface of a water container. The surface tension at the contact line near the rear of the toothpick 
is smaller than that near the front, leading to forward propulsion. The toothpick length is 6.5 cm. 



has been proposed [40, |41| • The combination of a chemical and hydrodynamic field has further been found to lead 
to interesting coupling between different bodies 421. and between bodies and boundaries [i^. Similar physics 
dictates the motion of camphanic acid boats |44l. |4£ , phenanthrolinc disks [i^, , ethanol-soakcd gels on water 
surface [13 and mercury drops near oxidant crystals 49| . 

In this paper we consider what is perhaps the simplest model of Marangoni propulsion, namely the case of a thin 
circular rigid disk located at the surface of a viscous fluid f j|ll| . We assume that a portion of the perimeter of the disk 
has a prescribed concentration of insoluble surfactant, to which the rest of the disk surface is impenetrable. Assuming 
a linearized equation of state between surface tension and the concentration of surfactant, we derive mathematically 
in the limit of low Capillary, Peclet and Reynolds numbers the surfactant concentration fleld ( §111 Ap . the 



fluid velocity and pressure fields ( mil B[) and the propulsion speed of the disk ( mil Cp . Our results, and the regime 
in which all our assumptions are expected to be valid, are summarized in ^IVI In particular, we show that neglecting 
the viscous stresses acting on the disk as a result of Marangoni flows over-predicts the propulsion speed by 50%. 



II. MODEL PROBLEM AND SETUP 



The model problem we consider in this paper is illustrated in flgurc [21 A solid disk of radius a lies in the horizontal 
surface of otherwise unbounded, viscous, incompressible fluid whose coefficient of viscosity is fi. Fluid motion and 
edgewise translation of the disk are induced by a non- uniform surface distribution of an insoluble surfactant, of 
concentration F per unit area. Specifically, the rear part of the disk (the arc described by |(/)| between a and tt, where 
(j) is the polar angle in cylindrical coordinates) is assumed to have a prescribed concentration of surfactant, described 
by /(|(/)|) in a dimensionless form. The front part of the disk {\(j>\ between and a) is impenetrable to the surfactant 
and thus the normal derivative of F is zero on this circular arc. The surfactant is assumed to be insoluble in the 
fiuid, and thus we ignore its transport to and from the bulk. For a water-air interface, this would be an appropriate 
approximation for surfactants such as, e.g., oleyl alcohol or hemicyanine [soj . and more generally organic compounds 
with long-chain hydrocarbon tails and polar hydrophilic heads j51| . 

Provided we consider this problem after the initial transient in fluid flow and surfactant transport, the governing 
steady-state equation for surfactant concentration at the free surface [s^] is given by 



W^F = V-(v,F), 



(1) 




FIG. 2: Setup for our model of viscous Marangoni propulsion. A solid disk of radius a is located at the surface of a semi-infinite 
fluid of viscosity /i. The front part of the disk (|0| between and ct in polar coordinates) is assumed to be impermeable 
to surfactants. The rear part of the disk (|0| between a and tt) is assumed to have a prescribed concentration of insoluble 
surfactant, described by f{\4>\) in a dimensionless form, leading to forward propulsion of the disk at speed U . The free surface is 
assumed to be undeformed by the motion of the disk, and remains flat (limit of small Capillary numbers). The Peclet number 
for the surfactant is small enough to neglect all convective terms in the surfactant transport equation, while the Reynolds 
number in the fluid is small enough for the fluid motion to be described by Stokes flow. 



where D is the diffusion coefficient and v,, is the surface component of the fluid velocity v. This latter satisfies the 
continuity equation, 

V • V = 0, (2) 

and is assumed to be sufflciently small not only to allow use of the Stokes equations, 

AiV^v = Vp, (3) 

where p is the dynamic pressure, but also to furnish a zero Peclet number approximation to the surfactant field, Fq, 
by neglecting the right hand side of ([T]). In addition, the fluid surface is assumed to remain flat. In other words, we 
focus in this paper on the limit where all three Reynolds, Peclet and Capillary numbers are asymptotically small; an 
a posteriori derivation of the conditions under which these limits are valid is offered in ^IVI The zero Peclet number 
surfactant field Fq, satisfies Laplace's equation, and is fully determined by its prescribed values at the disk ( ijlll Ap . 
The resulting Marangoni stresses induce the lowest approximation to v f ijlll Bp . which enables the translation speed 
of the disk to be determined f ^illl Cp . We use cylindrical polar coordinates (or, </), az), with the origin at the center of 
the disk and z measured vertically downwards (see figure [2]). In what follows, both r and z are thus dimensionless. 



III. ANALYSIS 



A. Surfactant concentration 



Being a solution of the two-dimensional Laplace's equation, the surfactant concentration, Fq, in the surface z ~ 0, 
exterior to the disk, is of the form 

oo 

Too? = Co + 2 ^ c„r"" cosn0, (r > 1). (4) 

n=l 

The dimensionless coefficients {c„;ri > 0} arc determined by the mixed boundary conditions, 

Loa" = /(I0I) at r - 1, a< |0| < tt, (5a) 

— l = Oatr=l, O<|0|<a, (5b) 
or 
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in which / is a prescribed, non-constant function that is symmetric with respect to the (j) = ±7r direction. AppUcation 
of (jSap and (j5b[) to ^ yields a mixed boundary value problem involving cosine series, namely. 



cq + 2 ^ c„ cos 710 — f{(l>) {a < (f) <tt) 

n=l 
oo 

nc„ cos 71,0 = (0 < (/i < a). 

n=l 

Following the method given by (§5.4.3), we set 

oo 

2 c„ COS 7 



Co 



1 

cos-i; 



h{t)dt 



-y/cOS (/) — cos t 



(0 < (/)< a), 



(6a) 
(6b) 

(7) 



which extends the range of (|6ap to the complementary interval by using an Abel transform in anticipation of the 
y/a — (j) behavior as — > a from below. This immediately yields the Fourier coefficients 



1 



TT 



1 

cos —V 
2 



h{t)dt 



: COS nv dv 



f{v) cosnv dv 



{n > 0). 



(8) 



ly V COS V — COS t 

When dH]) is substituted into the summations can be evaluated in closed form; but first we need to relate h{t) to 
A change in the order of integration in ^ yields 



hit)dt 



cos -^v COS nv dv 



f{v) cos nv dv 



{n > 0). 



(9) 



Iq \/cOS V — COS t 

When these coefficients are substituted into the integral of (|6b|) between and - the integration constant being 
trivially zero - the common summation can be evaluated and the double integral reduced to a single integral, giving 
the Abel integral equation, 

(">> 

iO<(f>< a). 



h{t) dt 



2 1 
= — sm-c 

V cos t — cos TT 2 ' 

Inversion then gives, after evaluating the t-integral. 



f{v) dv 



COS( 



cosv 



2 d r f{v)sinh, 



TT dt 7a -y/cOS t — COS V 

The first two coefficients are thus given, from ©, by 



:dv, 



(0 < i < a). 



Co 



Cl 



V2J0 



— / h{t)dt + - / f{v) dv = — 



y/2 r /(w)sinii 



TT j ^ ^COS OL — COS V 



zdv. 



If" 1 I r 

—= / h {t) cos^ -t dt -\ — / f{v)cosvdv 

V 2 Jo 2 TT 

(1 — COS a + 2 COS v)dv. 



1 



^ f^v)smjv 



(10) 



(11) 



(12) 



(13) 



nV2 J a V cos a — cos v 

More significantly, the solution ((TT|) enables the integral in ([5]) to be expressed in terms of the prescribed /(0) by 
evaluating the t-integral to give 



h{t) dt 



^f{v)sm^vdv Icosu— cos a 



(0 < u < a). 



(14) 



\/cos u — cos t TT COS u — cos v V COS a — cos v 

Now the substitution of ([8]) into (jlj and evaluation of the summations yields a closed form expression for the 
surfactant distribution, namely. 



1 



1 



cos —u 



+ 1 — 2r cos(0 + v) + 1 — 2r cos(0 — v) 

f{v)sm^v I cos u — cos a 



dv 



-dv 



cos u — cos V V cos a — cos w 



- 1 



1 



+ 1 — 2rcos(0 + u) + 1 — 2rcos(0 — u) 



du, 



(15) 
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valid for r > 1. Note that the case where f{(j)) = /oi a- constant, yields a uniform surfactant concentration around the 
disk edge, and hence no asymmetry to generate motion. 



B. Pressure and velocity field 

The Marangoni stresses induce a fluid motion via the boundary condition for tangential stresses 

oz 



(16) 



where 7 is the surface tension, is the gas constant, Ta is the absolute temperature, and the linearized equation of 
state, 7 = 70 - 7^TA^ [H, [Hi, is invoked. 

The pressure p has the Fourier mode expansion. 



p = Po (J^, -z) + 2 ^ P„ (r, z) cos n(^, 



(17) 



in which, since ^ and ([3]) imply that p is Laplacian, the coefficients have Hankel transforms of type 

2[i 



P„ = / kA„{k)e^''''J„{kr)dk, {n > 0). 
a Jo 

It is then advantageous to write the velocity field in the form 

00 

V = er.S'o(r, z) + y~^[er cos n(j> + smn(j>]Sn{r, z) 

n=l 

00 00 

+ ^[e^ cos n(j) — smTi(j)]Qn{r, z) + ez[Wo{r, z) + 2 ^ W„(r, z) cos ncj)]. 

n—l 71=1 

This ensures that the modal structure is preserved on substitution, with PT)) . into ([3]) which gives 
92 Id 1 92 



E 

71=1 

E 



Gr COS nm + sm nd 



Gr cos nm — sm nc 



a 

00 







" 92 


1 9 


9r2 


r dr 


" 92 


1 9 


9r2 


r 9r 




00 



92 ■ 

9z2 


S71 


92 ■ 


Qn 


9z2 




n2 


92 




9z2 



COS no 



, , ^e^ COS n0 + erf, sm n0J — P„ 

or -"^ — ' \ or r 

n— 1 ^ 



/ dP n \ I dP ^ Qp 1 

cos 710 - Grf, sin 710] I — ^ H — Pn I + Wo + 2 — ^ cosri^ 

n=l ^ ^ \ n=l J 



(18) 



(19) 



(20) 



With Pn given by ([T8| , it readily follows that the three sets of partial differential equations in ([20| have solutions 
of the form 



[Dn{k) + kzAn{k)]e-''\Jn+i{kr)dk, {n > 0), 
[Bn{k) - kzAn{k)]e~'''Jn^i{kr)dk, {n > 1), 



Wn ^ I kzAn{k)e~''''Jn{kr)dk, {n > 0), 
'0 



(21a) 
(21b) 
(21c) 
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in which the complementary function in (|21c|) is chosen according to the assumption of no normal flow in the plane 
{z = 0) of the disk and surfactant. Substitution of into ([2]) gives 

i d dWo ^ 



1 d 



^-1 d ( Qn 



dz 



(n> 1), 



whence 



Do(fc) + ^o(A:) -0, I?„(A:)-B„(A:) + 2A„(fc) =0, (fc>0,n>l). 



(22) 



(23) 



These functions are determined by imposing conditions on the tangential velocity components in the plane z = 0. 
The uniform translation of the disk, due to the Marangoni stresses is expressed by 



V = Ugx at z = 0, 7' < 1, 
where U is to be determined by a subsequent force balance (see ^III C|) . When substituted into ([1 

Qi^U, Qn = 0{n>l), = (n > 0) at z = 0,r < 1, 



(24) 

implies that 
(25) 



As in [Ha, [531, the Besscl functions in pia^ and (j21b|) arc changed to J„ by integration with respect to r, aided by 
the classic recurrence relations [ssl ] . The arbitrary constants that multiply a non- negative power of r are retained and 
chosen below to eliminate square root singularities in the velocities. Thus our required forms of the integral conditions 
on the interval r < 1 are 



Dn{k)k-^ Jn{kr)dk = d^r" , (r < l,n > 0), 
B„{k)k~^J„{kr)dk = ^rSi^, (r < 1,77. > 1), 



(26a) 
(26b) 



where (5,„„ denotes the Kronecker delta. 

Integral conditions on the interval r > 1 are obtained by applying the Marangoni stress condition, (jl6|) . Equation 
(|4]) shows that only the functions {S'„;77, > 0} in (|T9)) arc directly forced by the surface tension variations in (|16p . 
Replacing J„±i by J„ as above, our required forms of the integral conditions become 



[Do{k) ~ Aoik)]Jo{kr)dk = 0, (r > 1), 



[Ai(fc)- A„(fc)]J„(fcr)dfc = 



TZTa 
2"' 



(r>l,77>l). 



[Bn{k) + An{k)]Jn{kr)dk = br,r-'' (r>l,?i>l). 
We now use (|23p to eliminate {A„; tt, > 0} from p7a|) - ()27c|) and obtain 



Doik)Jo{kr)dk = 0, (r > 1), 



Dn{k).Jn{kr)dk = 
B„{k)J„{kr)dk = 



(r > l,7i> 1), 
(r > 1, 77 > 1). 



(27a) 
(27b) 
(27c) 

(28a) 
(28b) 
(28c) 



Equation (|26ap with (|28ap . (|28bp and (|26bl) with (|28cp define two infinite sets of mixed boundary value problems of 
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a type discussed by [s^ (§4.3). The continuation of the integrals in (|26a|) and (j26bp is given by 

2do r^/'' du 



Do{k)k-\Jo{kr)dk 
Dn{k)k^^ Jn{kr)dk 

Bn{k)k-^Jn{kr)dk 



(r > 1), 



r(i/2)r(n + i/2) Jo \/r^ 



^^^^ 



r(n- l/2)(r2 - 1)1/2 



4r(l/2)r(n)r" 



Jin 

TT 



35„ 



r(n- l/2)(r2 - 1)1/2 
4r(l/2)r(n)r" 



(r > l,n> 1), 



(r > l,ri > 1). 



(29a) 
(29b) 

(29c) 



The arbitrary constants are finally determined by eliminating edge singularities from the surface velocity compo- 
nents. From (|21ap and (|29ap . we have 



^o(r,0) 



Do{k)Ji{kr)dk 



d 



2do d 



— — arcsin - = 0, 

TT dr 



Da{k)k-^ Ja{kr)dk 
(r > 1), 



(30) 



provided that do = 0. 
From (PTa|) and ((29b| . 



f°° d [ \ f°° 

5„(r,0) = / Dnik)Jn+i{kr)dk^-r"— — D,-,{k)k-^ J„{kr)dk 
Jo dr [r" Jo 



2d„r{n + 1) 



„-(«+!) 



r(l/2)r(n+l/2) (r2- 1)1/2 



r"r(n- 1/2) d 
4r(l/2)r(n) d? 



(,.2 _ 1)1/2 



r(n + l/2) (r2- 1)1/2 
2r(l/2)r(n) ^^^^+1 



(r>l,n>l), 



provided that 



From ((2Tb| and (P^c)) . 



2d„r(n + 1) 
r(l/2)r(n+l/2) 



3 ^c„ + bn 



r(n- 1/2) 
4r(l/2)r(n)' 



(n>l). 



(31) 



(32) 



gn(r,0) 



/ Bn{k)Jn-i{kr)dk = ——r'' Bn{k)k-^Jn 
Jo ^ dr I Jq 

(5i„ arcsin [—) , (r>l,n>l), 



(kr)dk 



(33) 



provided that 



5-c„ + 3o„ = — din, ['n>l). 



(34) 



On using ([32]) and (|34p to eliminate the coefficients {d„,6„;n > 1}, the surface velocity is given from ([T9| by 

. ^, (r2- 1)1/2 



2f/ 

TT 

3//a2 



(e^ cos ( 



cos nm + Bj, sm na 



3r2 

r(n + l/2) (r2- 1)1/2 
r(l/2)r(n) r"+i 



COS + ga sm < 



(r>l). 



(35) 
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Note that = Ue^ + ©[(r^ - l)i/2] 

as r — >■ 1+. giving the expected continuous velocities. The summation in (|35p 

can be expressed as 

(r^ - l)i/2v 

and then, by reference to (1 — x)~^^^ , evaluated as in ([15]). 



n=l 



, r(n + l/2) 

'r(i/2)r(n + i) 



(r > 1), 



C. Marangoni propulsion speed 



The disk speed U is found by enforcing that the viscous drag force, 

2TTfie.^ [ ^(r,0)rdr = -27r^e^ / \Bi(k) + Ai(k)]Ji(k)dk, 
Jo oz Jo 

is balanced by the force on the disk due to surface tension, 

{j)r=i cos (j)d(j) = -2iTeJ^^ci, 



(36) 



(37) 



in which (|19p . (j21bp . (fT6)) . (|4]) have been used (recall that ci is given explicitly bv[T3|. Thus U is determined by the 
equality 



[Bi{k) + Ai{k)]Ji{k)dk ^ ^ci. 

The continuation of the integral in (|27cl) is given, for n = 1, by 

(3C/-2di)r bi 



[Bi{k) + Ai{k)]Ji{kr)dk = 



7r(l-r2)i/2 



1 - 



1 

3^ 



8U TZTa 



-ci 



(1 -r2)i/2 



(38) 



(39) 



(r<l). 



TT /iO^ 

and substitution into psp allows us to obtain the final expression for the Marangoni propulsion speed as 

For the particular case, /(0) = — cos0, we get, on substitution in p^ . (|13p . 

Co = -(1 — cosa), ci = — -(1 + cosa)^. 
2 8 



(40) 



(41) 



The positive surfactant concentration at the rear of the disk induces a lower surface tension there and hence a positive 
speed U. Note the limit values, Co(0) = and Ci(0) = —1/2, corresponding to Ta^ = —r~^cos4> in Also 
co{tt) = 1 and ci(7r) = 0, corresponding to a uniform surfactant of magnitude /(tt), which is the a — > tt limit of 
([T2|). Similarly, the particular choice /(</>) = cos2(/) yields ci = — 1(1 + cosa) sin^ a, while /(0) = — cos3(/) yields 
Ci = — ^(1 + cosa)(l — 5 cos a) sin^ a. 



IV. DISCUSSION 



In this paper we consider a simple model setup for Marangoni propulsion: a thin rigid circular disk located at the 
fiat free surface of a viscous fiuid and releasing insoluble surfactants with a prescribed concentration along its edge. As 
seen above, although both the geometry and setup are quite simple, the mathematical details of the derivation can be 
quite involved, and point to the richness of the physical problem. Physically, in order for the disc to move with finite 
speed, Eq. (|40p shows that the surfactant concentration cannot be uniform. From a thermodynamic point of view, the 
disk motion is thus an intrinsically non-equilibrium effect, where a non-uniform surfactant chemical potential needs 
to be maintained along the disk edge to drive the Marangoni propulsion. 
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The steady-state analysis in the paper makes three distinct physical assumptions, namely the limit of small Capillary, 
Peclet and Reynolds numbers. The result for the propulsion speed, PO)) . allows us to derive the appropriate velocity 
scale for the fluid, namely U ~ TZTaci/^o^. With a natural length scale L ^ a, wc can then use this result to derive 
the regime under which our three assumptions are expected to be correct, namely 

fiU UTaci 

Ca = = — < 1, (42a) 

70 loa^ 

An interesting observation from our theoretical results concerns the role of Marangoni stresses on the disk motion. 
An approximate method to derive (j40p would consist in neglecting the stresses on the disk arising from Marangoni 
flow. In that case, U would be obtained by balancing the direct forcing from surface tension at the contact line, 
F = —2TrTZTACi/a ((37l with ci from [T3)) by the viscous drag on the bottom half of the disk, F = 16fiaU/3 \E^ . 
leading to propulsion at the speed U = —STrTZTACi/Sfxa^ . This result, although it has the correct sign and order of 
magnitude, actually over-predicts the propulsion speed (j40[) by 50% and thus both direct surface tension forcing at 
the contact line and Marangoni flow play important role in determining the steady-state Marangoni propulsion speed. 
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